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' Abstract 

I The gauge-invariant three-gluon vertex obtained from the pinch technique 

0^ ■ is characterized by thirteen nonzero form factors, which are given in complete 

, generahty for unbroken gauge theory at one loop. The results are given in d 

I dimensions using both dimensional regularization and dimensional reduction, 

' including the effects of massless gluons and arbitrary representations of mas- 

sive gauge bosons, fermions, and scalars. We find interesting relations between 
I the functional forms of the contributions from gauge bosons, fermions, and 

I ' scalars. These relations hold only for the gauge-invariant pinch technique ver- 

tex and are d-dimensional incarnations of supersymmetric nonrenormalization 
r~| ! theorems which include finite terms. The form factors are shown to simplify 

for N = 1,2, and 4 supersymmetry in various dimensions. In four-dimensional 
non-supersymmetric theories, eight of the form factors have the same functional 
form for massless gluons, quarks, and scalars, when written in a physically motl- 
ey I vated tensor basis. For QCD, these include the tree-level tensor structure which 
has prefactor /3o = {llNc — 2Nf)/3, another tensor with prefactor 4Nc — Nf, and 
six tensors with Nc — Nf . In perturbative calculations our results lead naturally 
to an effective coupling for the three-gluon vertex, a{kf, /c|, A:|), which depends 
on three momenta and gives rise to an effective scale Q'^^j!(ki,k2,k'^) which 
governs the behavior of the vertex. The effects of nonzero internal masses M 
are important and have a complicated threshold and pseudo-threshold struc- 
ture. A three-scale effective number of flavors A''ir(A;^/M^, fc^/M^, fc|/M^) is 
defined. The results of this paper are an important part of a gauge-invariant 
dressed skeleton expansion and a related multi-scale analytic renormalization 
scheme. In this approach the scale ambiguity problem is resolved since physical 
kinematic invariants determine the arguments of the couplings. 



^Work supported by the Department of Energy under contract number DC-AC02-76SF00515 



1 Introduction: Gauge-Invariant Green's Functions 



The main purpose of this paper is to analyze the structure of the gauge-invariant 
three-gluon vertex calculate the fourteen form factors at one loop, and outline 
some of the phenomenological applications. Before proceeding, it is worthwhile to 
review the motivation and current status of gauge-invariant Green's functions. 

In the conventional formulation of gauge field theories, the manifest gauge-invariance 
of the original action is lost upon quantization, simply because one has to fix a gauge 
in order to perform calculations. Generically, Greens's functions are gauge-dependent 
and thus not physical by themselves. Only the particular combinations of Green's 
functions which form physical observables must be gauge- invariant. In many the- 
oretical studies, however, one would like to consider individual Green's functions 
and extract physical meaning from them For example, studies of the infrared 
behavior of gauge theory using Dyson- Schwinger equations [3] often rely on gauge- 
dependent truncation schemes which one hopes are not too brutal. The existence of 
gauge-invariant two-point functions is crucial for defining meaningful resummed prop- 
agators jlj, particularly near threshold, for the construction of effective charges [Sj, 
for a postulated dressed-skeleton expansion of QCD [B], and for justifying renormalon 
analyses jTj. 

Thus there is strong motivation for gauge-invariant Green's functions with phys- 
ical content. We will now briefly discuss the relationship between three different 
approaches to gauge-invariant Green's functions : (1) the Pinch Technique (PT), (2) 
the Background Field Method (BFM), and (3) the -k effective Lagrangian scheme of 
Kennedy and Lynn. All three approaches will lead to the same Green's functions. 

The pinch technique (PT) was first constructed by Cornwall [2j in order to study 
gauge-invariant Dyson-Schwinger equations and dynamical gluon mass generation, 
but the approach is much more generally applicable. In the PT approach, unique 
gauge-invariant Green's functions are constructed by explicitly rearranging Feynman 
diagrams using elementary Ward identities (WI) as the guiding principle. Longitudi- 
nal momenta from triple-gauge-boson vertices and gauge propagators inside of loops 
hit other vertices and thus generate inverse propagators (via WI's), which, in turn, 
cancel (or pinch) some internal propagators. In this way, certain parts of Green's 
functions are reduced to parts of lower ra-point functions, and should properly be 
included in the latter. 

As an example of the PT, consider the gluon (or massive gauge boson) self-energy. 
The conventional self-energy is gauge- dependent and physically meaningless by itself. 
However, when embedded in any physical process, there will be associated parts of 
vertex and box graphs which undergo the reduction described above and thus have 
the same tensor and kinematic structure as the gluon propagator. These pinched 
parts are then added to the conventional gauge-dependent self-energy, yielding a 
gauge-invariant self-energy and gluon propagator that has the correct asymptotic 
UV behavior dictated by the renormalization group equation. The resulting two- 
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point function has numerous positive attributes |1] [HI [H] [ID] , including uniqueness, 
resummability, analyticity, unitarity, and a natural relation to optical theorem, from 
which it can also be derived ji] ^l] . 

Resumming these two-point functions leads to physical effective charges, ala 
Grunberg|12j. which can be extended to the supersymmetric case and leads to an 
analytic improvement of gauge coupling unification with smooth threshold behavior 

This method has been applied to a variety of Green's functions [THjpTj . 
with applications to electroweak phenomenology [IHl[ini- In particular, the gauge- 
invariant three-gluon vertex was first constructed in P] to one-loop order, where the 
authors showed that the vertex satisfies a relatively simple abelian-like Ward iden- 
tity. However, the integrals were not evaluated, so that little could be said about the 
individual form factors except that the UV divergent term in the tree level tensor 
structure is correct. The main motivation of this paper is to extend this work by 
evaluating the integrals for the fourteen form factors, and expressing the results in a 
convenient tensor basis for phenomenological applications. In doing so, an interesting 
structure emerges, in which the contributions of gluons(G), quarks(Q), and scalars(S) 
are intimately related. These relations are closely linked to supersymmetry and con- 
formal symmetry, and in particular the iV = 4 non-renormalization theorems. For all 
form factors F in dimensions d, we find that 

FG + AFQ + {lQ-d)Fs = Q, (1) 

which encodes the vanishing contribution of the = 4 supermultiplet in four dimen- 
sions. Similar relations have been found in the context of supersymmetric scattering 
amplitudes [201 In Appendix E, the effects of internal masses are discussed, and 
the above sum rule becomes modified 

Fmg + 4Fmq + (9 - d)FMs = 0, (2) 

for internal massive gauge bosons (MG), fermions (MQ), and scalar (MS). The exter- 
nal gluons remain massless and unbroken, so the internal gauge bosons might be the 
heavy X^,y^ bosons of SU{5), for example. In |221; supersymmetric relations were 
found for electroweak gauge boson four-point scattering amplitudes. 

The PT method has been explicitly extended beyond one- loop j^^f^f^, has 
recently been proven to exist to all orders in perturbation theory [211 Elij li28j |^29jj , and 
interestingly, each Green's function is equal to the corresponding Green's function 
of the Background Field Method (BFM) in quantum Feynman gauge = 1, a 
result suggested in [20] [211 • Heuristically, this is due to the fact that there are no 
longitudinal (pinching) momenta in the gauge propagator or the elementary vertices 
in this special gauge. 

The Background Field Method (BFM) [221 constructs manifestly gauge invariant 
Green's functions in the following way. First, the field variable(A) in the path integral 
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is separated into a background(B) and quantum((5) field, A = B + Q. Only the 
quantum field Q propagates in loops, since it is a variable of functional integration. In 
contrast, the background field B appears only in external legs. By judiciously choosing 
the gauge-fixing function, one arrives at an effective action which remains manifestly 
invariant under background field gauge-transformations SB"^ = —f°''"^u!^B'^ + ^d^u"'. 
Furthermore, derivatives of the BFM effective action with respect to the background 
field B yield the same IPI Green's functions as the conventional effective action with 
a nonstandard gauge-fixing. Thus, it can be shown that the correct S-matrix is 
obtained by sewing together trees composed of IPI Green's functions of B fields. In 
doing so, one can fix the gauge of B, which propagates only at tree level, independently 
of the gauge fixing of Q. For example, convenient non-covariant gauges might be used 
for the trees while BFM Feynman gauge = 1 (BFMFG) can be used for the loops. 

The correspondence between the FT and BFM is not surprising, since the BFM is a 
formulation of gauge theory where Green's functions of the gauge field are manifestly 
(background) gauge-invariant. Although this is true for all values of the quantum 
gauge-fixing parameter ^g, it is only for the special value = 1 that the BFM Green's 
functions also have the correct kinematic structure of the irreducible FT Green's 
functions. Alternatively, it has been shown that applying the FT algorithm to 
the BFM for 1 leads back to the canonical (^g = 1) FT Green's functions. 

Finally, in the * scheme of Kennedy and Lynn j^SI; a gauge-invariant effective 
Lagrangian was constructed for electroweak four-fermion processes by explicitly re- 
arranging the one loop corrections. As in the pinch technique, vertex parts must be 
added to would-be two point functions to yield genuine two-point functions. One 
particular motivation is that fact that the photon acquires a spurious mass from its 
mixing with Z°, n^^(g^ = 0) 7^ 0, unless the correct vertex parts are added. The 
resulting effective charges, and s^(g^) are in fact precisely equal to the corre- 

sponding pinch-technique effective charges at one loop, including all finite terms and 
threshold dependence |S]. 

Thus, all three methods for constructing physical gauge-invariant Green's func- 
tions lead to the same results, which in the this paper will be referred to as either FT 
or FT/BFMFG Green's functions. 

The organization of this paper is as follows. In section 2, we will discuss the 
general structure of the gauge-invariant three-gluon vertex, which is constrained by 
the Ward identity and Bose symmetry. Two convenient tensor bases and their rela- 
tion are discussed. In section 3, the main results of this paper are given. First, the 
nontrivial supersymmetric relations between the gluon, quark, and scalar contribu- 
tions to each form factor are discussed. The explicit results for the form factors are 
given in two different bases for massless internal particles, with the full mass effects 
relegated to Appendix E. In section 4, we briefly discuss the phenomenological ap- 
plication to physical scattering processes, where we derive an effective coupling for 
the three-gluon vertex, 5(/cJ, and an effective scale, (^^^.^(/c^, /c|, both of 

which depend on three distinct gluon virtualities. In section 5, the phenomenolog- 
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ical effects of internal masses are discussed. A complicated threshold and pseudo- 
threshold structure emerges. Furthermore, a three-scale effective number of flavors 
A^i?(/c^/M^, /c|/M^, /c|/M^) is defined. Conclusions are given in section 6. In ap- 
pendix A, a brief outline of the calculational method is given, and some basic one- 
and two-point integrals are given. Appendix B is devoted to a thorough discussion of 
the massive triangle integral, and analytic continuations are given for each kinematic 
region. Appendix C collects some useful results for special functions. Appendix D 
explains the corrections to the form factors when a supersymmetric regularization is 
used. Finally Appendix E gives explicitly the corrections to the form factors arising 
from internal massive gauge bosons, fermions, and scalars. 



2 General Structure of the Three-Gluon Vertex 
2.1 Symmetries 

One of the most important aspects of the gauge-invariant three-gluon vertex discussed 
in this paper is the relatively simple Ward identity it satisfies, which has the same 
form as the Ward ID satisfied by the tree level vertex. This was proven at one-loop 
in the original paper by Cornwall and Papavassiliou pDj using the explicit one-loop 
result, which is the gluon part of Eqs. ()17|18|l below. It is straightforward to show that 
the fermion and scalar parts also satisfy the same Ward identity (just as in QED). 
Furthermore, the equivalence of the BFMFG and PT to all orders j^H] allows one to 
write the Ward identity satisfied by the three-gluon vertex to all orders as 

(3) 

plus two other equations which are cyclic permutations. The transverse tensor t^^^p) = 
P^Qiiu — PfiPu comes from the tree level term. Here all momenta are defined to be in- 
coming and all labels are defined in counter-clockwise fashion, as shown in Fig. 1. 
This Ward identity represents a great simplification compared to the usual Slavnov- 
Taylor identities satisfied by the conventional gauge-dependent three-gluon vertex, 
which involves the gluon propagator, the ghost propagator, and the ghost-ghost- 
gluon vertex function. The self-energy function in the above equation is not the usual 
gauge-dependent self-energy, but rather the gauge-invariant pinch technique self en- 
ergy, which is the only self energy discussed in this paper. An immediate consequence 
is that the longitudinal (L) part of the vertex, defined as the part which contributes 
to the above Ward ID, must have only the antisymmetric color factor f"''"^ so long as 
gluons conserve the color charge, 11^^ (g) = 5'^*n^j^(q'). As far as we know, the trans- 
verse (T) part of the vertex (defined by P3^Tl^^^2^^{pi-, P2-, Ps) = 0) is not required to 
be proportional to f"-'"^, but may in principle contain ci"'"^ terms. Nevertheless, no 
such terms appear at one or two loop order, and so in the subsequent discussion we 
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Figure 1: The notation and loop momentum routing used throughout this paper. The 
internal particle could be a gauge boson, ghost, quark, or scalar. 



take 

Cm.M3(P1'P2,P3) = r"'T,,,,,M,P2,P3), (4) 

in terms of which the Ward identity becomes 

Bose symmetry, the fact that 3 identical particles are entering the vertex, and the 
properties of f"'^'^ imply definite properties of T^j^^^^^{pi,p2,P3) under the interchange 
of labels. In particular, defining the five elements of the permutation group S3 to act 
by 

f (/^2,P2) (/i3,P3) \ n - f if^3,P3) ^ if^uPl) \ 

V (/^3,P3) - (/i2,P2) J V (/^I'Pl) ^ (/^3,P3) J 

and 5(321 = fi-fag one finds that (5-123, 5'32i, 5'i2, 5'23, 5'3i) yields (+,+,-,-,-) when 
acting on T^^f,^^^{pi,p2,P3). 

The nonabelian nature of the permutation group 1S3 prevents one from finding 
a basis in which all of the tensors are eigenstates of all of these operators. Thus, 
aesthetic and physical principles must guide us in choosing convenient bases, two of 
which are discussed momentarily. 
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2.2 Two Convenient Bases 



Let us consider the most general tensor structure. The three index Lorentz covariant 
tensor must be constructed out of the metric g^i, and the momenta (pf). Since the 
momenta are not independent, pi + P2 + Ps = 0, simple combinatorics implies that 
there are in general 14 independent tensor components, 6 of which have one power 
of momenta and also the metric, and 8 of which have 3 powers of momenta. Many 
different basis choices can be made, although we will use essentially two. 

In the subsequent discussion, some efficient notation will prove useful. This is 
summarized in Table 1. 



Table I. Definition of tensor abbreviations 



001 = 5'a.im2Pim3 


002 = ^^i;,2P2m3 


200 = ^^2^3P2mi 


300 = g^,2t,3P3^,l 


030 = fi'^3^iP3A.2 


010 = fi-^a^^iPi^^ 


211 = P2^,,Pl,^2Pl^^3 


212 = P2^,^Pl^,2P2^,^ 


232 = P2^,,P3t^2P2t^3 


332 = P3a.iP3m2P2/.3 


331 = P3/.iP3m2Pim3 


311 = P3a.iP1m2P1M3 


312 =P3^iPi^2P2m3 


231 = P2,,iP3t,2Plt^3 



Thus, each tensor is rewritten as a 3 slot object, where slots correspond to 
/xi,/X2,/i3 in that order, and the content of each slot is either '1','2', or '3' to rep- 
resent momentum Pi,P2,P3, or a '0', which must occur in pairs and represents that 
those two indices are connected by the metric tensor. 

The most naive thing to do would be to just eliminate one momenta, say ps = 
-P1-P2 and use the following 14 basis tensors : 100, 200, 010, 020, 001, 002, 111, 112, 
121, 211, 122, 212, 221, 222. This is not very useful since the explicit Bose symmetry 
between the three gluons has been broken, and thus delicate relations between the 
form factors will have to enforce it. 

The ± Basis 

A more natural choice is obtained by starting from a manifestly symmetric, but 
redundant basis, which has 36 possible basis tensors, 9 with one power of momenta, 
and 27 with three powers of momenta. As a step towards our final basis, we find it 
convenient to eliminate all such tensors with momenta Pi^ ,P2^, or pg^, i.e. anything 
with 1 in the first slot, 2 in second slot or 3 in the third slot. This yields the 14 basis 
tensors 001, 002, 200, 300, 030, 010, 211, 212, 232, 332, 331, 311, 231, 312, which 
are shown in Table 1. Note that under the action of (7123 we have 200 — 030, 300 — >■ 
010,211 232,311 212, etc. Also, notice that 200 and 300 are interchanged 
by the action of (723, while 211 and 212 are interchanged by the action of gu, etc. 
Thus, it is convenient to take appropriate linear combinations such that one of these 
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interchange operators is diagonal for each tensor. Such basis tensors are 



ai2 = (00-) = 001 - 002, 023 = (-00) = 200 - 300, Ogi = (0-0) = 030 - 010 

6^2 = (00+) = 001 + 002, 623 = (+00) = 200 + 300, 631 = (0 + 0) = 030 + 010 

C12 = (+ + -), C23 = (- + +), C31 = (+-+) 

di2 = ( +), d23 = (+ 4i = (- + -) (7) 

h = ( ), s = (+ + +), 

where the notation means (±±±) = (2 ± 3, 3 ± 1, 1 ± 2), so that (+H — ) = 231 — 
232 + 211 — 212 + 331 — 332 + 311 — 312, etc. The subscripts are chosen because, for 
example, is an eigenstate of gi2, etc. 

Suppressing indices and momentum dependence, the three-gluon vertex is then 
written as 

r = {Auciu + i?i2&i2 + C12C12 + Di2di2 + perms) + 55 + Hh, (8) 

where the lower case letters represent the basis tensors, while the upper case letters 
are the form factors, which depend on p\,p% and p\. In addition to indicating which 
basis tensors they are associated with, the subscripts on form factors also indicate 
the ordering of momenta in the arguments. For example, A12 = ^(^1,^21^3)5 ^23 = 
A{p2,pI\p'1), A31 = A{pI,p1\pI), and the first two arguments are either symmetric or 
antisymmetric. The behavior of these form factors under 1S3 can be inferred from 
the behavior of the basis tensors under 1S3 (which will be discussed momentarily), 
along with the overall requirement for the vertex given below Eq. El One finds that 
A{x,y\z) = +A{y,x\z), and thus ^412 = A21, etc. Similarly, B12 = — -B21, Cu = C21, 
and D12 = —D21. H is totally invariant under the interchange or permutation of any 
momenta, while S goes to —5* under any interchange of momenta, but is invariant 
under a cyclic permutation (7123. 

It is straightforward to see that under the action of the permutation operator 
(^7123) these fourteen basis tensors are organized into four triplets, {012, ^23, 031}, 
{^12,^23,^31}, {ci2, C23, C31}, {(ii2, (^23, c?3i}, as well as h and s. The latter two are 
eigenstates of all five operators. 

Consider the properties of {012,023,031} under the permutation group. It is easy 
to see that under the action of any element gi, we have 

ai2 \ / ai2 \ 

023 ^ 9i\ «23 , (9) 
031 / \ ^31 / 
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with the matrices given by 



9l23 



9l2 




9321 



923 




1 

1 

1 



931 




(10) 



The transformation rules are identical for {cij}, and similar for {%} and {dij} 
with the only change being that there is no minus sign in the three interchange 
operators 5-12, 6-23 and 5(31. 

The basis constructed above (Eq.©) will be called the ± basis. As discussed 
later, this basis is the most convenient for phenomenology and furthermore the form 
factors exhibit particularly simple relations between the gluon, quark, and scalar 
contributions fEqs. ()45l55|) ). 

However, the ± basis as it stands does not contain the tree level tensor struc- 
ture. Thus, one is naturally led to diagonalizing the permutation operator (7123.^ 
Clearly, this is the most symmetric choice and, more importantly, one of the resulting 
eigenvectors is the tree level tensor structure. 

In the triplet representation of 1S3, (7123 is diagonalized by the similarity transfor- 
mation 




S 9123S 



9123 




:iii 



where A 



exp(^) 



:V3 



basis tensors and form factors 



A* are cube roots of unity. This results in new 





(12) 



This procedure is repeated identically for the {b,B), (c, C), and {d,D) basis tensors 
and form factors. 

Notice that do = 5-^1^2(^1 - ^2)^3 + 5-^2^3 - Ps)/.! + fl-Ms^ife - Pi)m2 is the tree 
level tensor, which is why the extra factors of were included above. 

^One can readily check that the only two operators in 53 which commute are (7123 and its inverse 
ff32i- Thus one can diagonalize these two, OR one of the interchange operators 512, <723, ffsi- 
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The transformation properties for the basis tensors cti and q are deduced from 
6^123 given above, 5^321 = diag(l, A, A), and 

/ioo\ /ioo\ 

5^12 = - 1 to = - A 5^31 = - A , (13) 
\010/ \0A0/ \0A0/ 

while the transformation properties of the form factors are deduced by demanding 
that behavior given below Eq.® is respected. For example, since a_ — Aa+ under 
g23, "we find A_ XA^ so that A_a^ — A_a_. For {b,B) and {d,D), the only 
change in the above is that there is not a minus sign in gi2,g23, and g^i. 

We have not touched h and s, as these are already eigenstates of all five operators 
{gi23, 9321, gi2, 923, 5'3i), with eigcuvalues (+ H ) and (+ + + + +), respectively. 

We will call the above constructed basis the symmetric ± basis. Note that any 
basis can be symmetrized in the same manner by diagonalizing the permutation 
operator (7123. The basis we started with in Eq.Q is motivated by (a) its simple 
and symmetric construction from only metrics and (pj+i±pj_i)^., (b) it is the most 
convenient basis for perturbative calculations, as will be discussed in section 4, and 
(c) The individual form factors have a relatively simple form, as will be discussed in 
section 3. 

The LT Basis 

For some theoretical studies, another convenient basis is determined by the dis- 
tinction between transverse (T) and longitudinal (L) tensors [SI] [SHI- The L tensors 
contribute to the Ward ID (or the more complicated Slavnov- Taylor ID for the gauge 
dependent vertex) while the T tensors satisfy homogeneous equations P3^r|i^)t2/i3 = 0. 
This is a very convenient basis for evaluating the loop corrections to the vertex, since 
the L and T parts separate, as described in Appendix A. 

The ± basis and the LT basis are complementary in the following sense. The 
± basis is constructed out of combinations of longitudinal (+) and transverse (— ) 
momenta, so that for example (00+) = fi'^i^alPi^s +^2/^3) = ~5'AtiM2P3M3 vanishes if the 
//3 index is contracted into a conserved current. Meanwhile, the LT basis distinguishes 
between parts of the vertex that do (T) and do not (L) vanish when dotted with 
longitudinal momenta. These straightforward relations to current conservation and 
Ward identities are essentially the reason these two bases are the most convenient to 
work with. 

In our notation, the vertex can be written in the LT basis as F = F^ + T^, where 

Tl = (v4i2ai2 + -B12&12 + C12C12 + perms) + Ss 

Tt = (Fi27i2 + perms) + :H7^, (14) 

and the bar distinguishes this LT basis from the ± basis defined above in Eq.((7)). 



10 



The relation between basis tensors is given by 



= 001 - 002 = ai2 (15) 
bi2 = 001 + 002 = 612 

C12 = 211-212-(pi.p2)(001 - 002) = ^(ci2-43 + 4i-/i)-(PrP2)ai2 

7l2 = (P1-P2)((P2-P3)001 - b3-Pl)002) - ((P2-P3)211 - (p3-Pl)212) 
{Pl-P2) 



2 



12 



+ ^ (^3(^12 - ^23 + 4l - /i) + (P? - P2) idl2 - C23 + C31 - S)J 



h = 231-312- (^(pi -Pa) (030 - 300) + (|)2-P3) (001 -010) + (p3 -Pi) (200 - 002) 

= ^{h + ci2 + C23 + C31) + ^ (^^3012 - (p? - ^2)^12 + Pld23 - {pI - pl)h23 

+ Pla^i - {pI - pDhi) 

s = 231 + 312 = ^(s + rfi2 + 43 + 4i), 

and we used {pvP2) = (Ps ~Pi ~P2)/2- This imphes the relation between form factors 

2 

A12 = A12- {pi-P2)Ci2-j[{Pi-P2)Fi2-H) 

2 2 

B12 = B^2 + ^^^({prP2)Fi2-H) (16) 



ri — ^ flT j_ ,Pi Ps-p ,P2 Ps-p 

^12 — ^ ^-H + <-'12 + — -T 12 H ^ 31 -\ 2 23 



2 / p^ p^ 

D12 = ^ (^5' + C23 — C31 + — F23 — —F3i-\ — F12 



TT _ ^ (tT 77 Ps-p 77 Pl-p 77 P2-p 
n — -^-H— O12 ^-ri2 — <-'23 ^-^23 — <-'31 2" 

S = lf5 + ^Fi2 + ^F23 + ^F3i 

4 V 2 2 2 

The unwritten form factors (A23, etc.) and basis tensors can be obtained trivially 
from the above equations by cyclic permutation ((7123) • In doing so it is useful to 
keep in mind the properties described under Eq.®, along with Fij = Fji, Aij = Aji, 
Bij = —Bji, Cij = Cji, while H and S have the same transformation properties as 
H and S, respectively. 
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3 Results for the Form Factors 



In this section, we will present the results for the form factors in arbitrary dimension 
d using dimensional regularization (DREG), for gluons in the adjoint representation, 
and massless quarks and scalars in arbitrary representations. The corrections due to 
supersymmetric regularization and massive fermions, scalars, and gauge bosons are 
given in detail in Appendices D and E, respectively. 

The gauge invariant vertex at one loop can be written as 



,2 



r2l2^3(pi'P2,P3) 

(17) 



\ f 

where the gluon (G), quark (Q), and scalar (S) integrals are 
d'^l 1 

+ 2(^2 + ^3)mi(^3 + h)fi2{h + h)fia ~ 8^l{9ti 

- 8^2(fi'M2M3P2/xi - fi'AtaMiPs/xg) " 8^3(9 ^l3^^lP3^l2 - 9^l3^^2P3^ll 



CmiM2M3 - J (^2vr)d/2/2]2 f r^Mi7*^^2,Pl, -h)'^^i,2aih, P2, "^1 )rf'^3/3 (^1 , PS, -^2) 



<5mim2M3 = J 7^^p^Tr[7^, |'37^2 /i7^3 /a] 



'S'AtiM2M3 - ~ /o„\rf ;2;2;2 (^2 + h)iii{k + ^1)^2(^1 + ^2)^3- (18) 



dH 1 
(2^/ff/i' 

The gluon contribution was first derived in P using the pinch technique(PT), and 
is equivalent to the vertex obtained in the Background Field Method in quantum 
Feynman gauge(BFMFG). The quark and scalar integrals come straightforwardly 
from the one loop triangle diagrams. The notation and routing of the integral are 
defined in Fig. ((H) such that /i = P2 + ^3, ^2 = P3 + ^i, h = Pi + h and the tree level 
vertex F^^^ and F-^ are defined as 

rSM2M3(Pl'P2,P3) = ^MiM2(Pl - P2)m3 + ^M2M3(P2 " Ps) + ^^3^ (PS " Pi), 



rJ),i^(/2,Pi, -^3) = '^Pip9^iii - 2pi75'Mi/3 - (^2 + h)i,^9M- 
All massless integrals can be reduced to two basic scalar integrals: 

dH 1 



J = j{p'1,pI,pI) = J 



{2nYllllll 



3 

dH 1 f dH 
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These functions, and the massive integrals which are considered later, are sum- 
marized in Appendices A and B, where J is written in terms of Clausen functions. 
In the following we will suppress the momentum arguments and write our results in 
terms of J, Ji, J2, Js- 

3.1 (Supersymmetric) relations between gluons, quarks, and 
scalars 

Before presenting the results for individual form factors, which are somewhat lengthy, 
we will discuss the relationship between the gluon(G), quark(Q), and scalar (S) con- 
tributions. We will end up finding relations similar to those found in the context of 
supersymmetric scattering amplitudes [20] [ZIl ISSj • For a generic form factor F, let 
us write the one-loop contribution as 

F = zg^ (CaFg + 2 ^ TfNjFg + 2 TsNsFs) , (21) 
/ 

where the coupling constant ig"^ and group theory factors have been pulled out. The 
standard notation is used, so that Ca = C*2(G) = Nc for SU{Nc), and Tr[tjtj] = 
TfS""^ . Thus, Fq stands for the contribution of one Dirac fermion in the fundamental 
representation of SU{Nc), or, due to a symmetry factor of | for Weyl fermions, the 
contribution of adjoint gluinos divided by N^. Similarly, Fs stands for the contribution 
of one complex scalar in the fundamental representation, or the contribution of a real 
scalar in the adjoint, divided by Nc- These identifications will be used shortly. 

After explicitly calculating the integrals in Eq.((THI), we noticed that the gluon, 
quark, and scalar contributions have a similar structure for each form factor. To 
make this explicit, define the following sums for form factor F : 

= - — - — ^-Fq + Fg 
EsGiF) = {d-2)Fs-FG. (22) 

Although the results for each form factor are often long, these sums are particularly 
simple, as can be seen in Eqs. (|42l45l55|) . For all form factors in any basis, it also turns 
out that 

{d-10)EsG = 8EQG (23) 

and Sqg is always proportional to d — 10. The above two equations and the results of 
Eqs. fl42l45l55|) can be used to determine the Q and S contributions to any form factor, 
given the gluon contributions written explicitly below. Furthermore, Eqs. ()22|23|l can 
be combined leading to 

Fg + 4Fq + (10 - d)Fs = 0. (24) 
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Considering the very different origins of each form factor fEqs. ()17ll8|) ). it is remark- 
able that they are related in such a simple manner. Note that no such analogous 
relation holds for the gauge dependent vertex |35] . 

This type of relation hints at supersymmetry. To further understand the content 
of these relations, we will consider various supersymmetries in d = 4. 

• N=l From the above definitions, it is clear that a vector superplet Vi (gluons 
plus gluinos) contributes ig'^Nc^Fc + Fq) = ig'^NcFy-^ to a generic form factor 
F, while chiral superplets contributes ig"^ Nf^ij^Fq + Fs) = ig^Ni^F^. By the 
sum rule Eg. ()24|) in = 4, we have Fy^ + 6-F<j, = 0. Thus any form factor can 
be written 

F = zg\N,Fv, + N^F^) = '^p\^=''^ Fy,, (25) 

where I3^~^^ = SN^ — |A^$ is the first /3 function coefficient. Hence the contri- 
butions of vector and chiral superplets have precisely the same functional form 
for each form factor. Furthermore, every form factor is proportional to /3o even 
though all but one of them are UV finite. 

• N=2 Here the vector superplet gives ig^N^Fc + 2Fq + 2Fs) = ig'^N^Fv^, 
while Nh hyperplets (a Weyl fermion of each helicity plus a doublet of complex 
scalars) yield ig'^Nh{FQ + 2Fs) = ig'^NhFh. The sum rule Eq.(E^ can be written 
as Fv2 + 2Fh = 0, and thus 

F = ig\KFy, + iV,F,) = =2)f^^, (26) 

where = 2iV, - Nh- 

• N=4 Here the vector superplet (the only multiplet allowed) contributes 2ig'^Nc{FG+ 
4Fq + QFs) = NcFy^, which is identically zero by the sum rule, which of course 

is a consequence of = 0. 

Thus, the similarities between form factors in (i = 4 are related to supersymmetric 
non-renormalization theorems. In particular, the exact conformal invariance of = 4 
implies that the gauge-invariant three-gluon Green's function is not renormalized at 
any order in perturbation theory. Furthermore, at one-loop order there are not even 
finite corrections, as refiected in Eq.(j2H). 

Analogous results hold for supersymmetry in ci 7^ 4. Here we must be careful, 
because in the sum rules and form factors presented in this paper we worked in d 
dimensions everywhere except in the traces over gamma matrices, where we used 
the conventional rule of dimensional regularization tr[7^7jy] = 4(yf^i,, and similarly 
for other traces. Properly working in integer valued d dimensions, we should use 
ti'b/xTy] = ds{d)g^^, where the spinor dimension of the gamma matrices is 

w ,^ f 2"'/2 \ f d even \ 

ds{d) = ^ 2(<^-i)/2 ) [d odd ) ■ (27) 
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Thus Fq^ 




Q is the contribution of a Dirac fermion in d dimensions, and we 



have 



Fg + 



16 



Fq^ + (10 - d)Fs = 



(28) 



ds{d) 



Rather than using Eq. (j28p . one can alternatively use Eq. (j24p and be sure to count 
fermion degrees of freedom in terms of d = 4 spinors. Thus the Weyl fermions of 
d = 6 and the Weyl-Majorana fermions oi d = 10 are composed of 2 and 4 Weyl 
fermions of four dimensions, respectively. From this, it is straightforward to show 
that d = 6, N = 2 and c? = 10, = 1 gauge theory give vanishing contribution to 
every form factor. For the d = Q, N = 1 case, one finds 



where jSo is determined from Eg. ()H7|) in d = 6. 

Note that it is not straightforward to analytically continue ds{d) into arbitrary 
non-integer d, which is the reason for the simple dimensional regularization rule. 
However, the sum rule expressed in Eq. ()24j) is an analytic function of d and thus rep- 
resents an analytic continuation of supersymmetric non-renormalization theorems to 
arbitrary d. This is intimately related to the existence of a supersymmetry preserving 
regulator, dimensional reduction (DRED), where vector degrees of freedom are kept 
in four dimensions while the integrals are still performed in d dimensions. Around 
four dimensions, d = A — 2e, we have Fq + 4Fq + (6 + 2e)Fs = in dimensional 
regularization, and we see that the e term plays the role of the so-called e-ghosts of 
DRED. We have calculated the form factors in DRED (see Appendix D) and verified 
that 



In the preceding discussion of supersymmetries in various dimensions we implicitly 
used DRED. 

The extension of these relations to the massive case is outlined in Appendix E, 
where the full effects of internal massive fermions (MQ), massive scalars (MS), and 
massive gauge bosons (MG) are included, and the sum rule becomes 



in DREG while in DRED the only change is 9 — (i is replaced by 5. Note that the 
external gluons remain massless and unbroken, so the internal massive gauge bosons 
might be the colored heavy gauge bosons arising in GUT models. The change of 
10 — (i in the massless case to 9 — (i in the massive case reflects the fact that massive 
gauge bosons "eat" one scalar degree of freedom. 

It should be emphasized that relations such as Eq. (j24p do not exist for the gauge- 
dependent three-gluon vertex [33], since the gluon contributions depend on the gauge- 
parameter, while the quarks and scalars do not. Indeed, it is uniquely the pinch 



F = -^(3oFv, (3o = 2N, - N^, 



(29) 



Fg + 4Fq + 6Fs = (DRED). 



(30) 



Fmg + 4Fmq + (9 - d)FMs = 



(31) 
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technique (or equivalently BFM in quantum Feynman gauge) Green's function which 
satisfies this homogeneous sum rule. For example, calculating in the BFM with 

7^ 1, leads to a nonzero RHS of Eq. (j24|l . 

Since the sum rule applies to all form factors, one finds 

G^^l^l2^^3 + ^Q^|.^^^2^^3 + (^^ ~ d)^^J■l^J.2^J.■i = 0' (32) 

which is remarkable given the expressions in Eq. ()18|) . One can explicitly show this by 
performing the trace in Eq. ljlSp and some tedious algebra to rearrange the T'^F'^r^ 
term. This can also be seen in the second order formalism of the BFM j^Hl |2I] [SH] • 

A similar relation holds for the one-loop gauge-invariant (pinch-technique) gluon 
two-point function in d dimensions, 

f 

(33) 

where from Eqs. ()3(ip37|l below we find 

G^..,2 + 4Q^iM2 + (10 - d)S^,^, = 0. (34) 

Unfortunately, analogous relations do not hold for higher gauge-invariant gluon 
n-point functions. This is essentially because the color and spacetime indices mix 
nontrivially. However, inhomogeneous relations of the form 

G + 4g + (10 -rf)^ = simple (35) 

still hold 1201123 EH]' where "simple" means an integral with fewer powers of loop 
momenta in the numerator. In the four-gluon case, this is just a simple scalar integral 
with no powers of loop momentum in the numerator. These loop momentum counting 
rules have been derived in the second order formalism, which is reviewed in jSE] • Note 
that the Ward ID for the four-gluon vertex J7j relates it to the three-gluon vertex, 
and thus the longitudinal parts of the four-gluon vertex must satisfy the homogeneous 
sum rules Fq + 4Fq + (10 - d)Fs = 0. 

It is interesting to see if extensions of these sum rules apply to two-loop calcula- 
tions, where the supersymmetric Yukawa vertices must be taken into account. As a 
first application, the two-loop pinch technique gluon self-energy has been calculated 
including finite terms. Interestingly, the finite terms do not vanish for = 4 SUSY, 
so it appears that the homogeneous sum rule in Eq. ()24j) does not have a counterpart 
at two loops. In any case, the finite parts of the two loop result allow for an im- 
proved extraction of the FT couplings from data as well as giving the three loop beta 
function. This calculation will be reported elsewhere [HTj . 

Now explicit expressions for the form factors will be given, first in the LT basis, 
and then in the ± basis. 
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3.2 The longitudinal form factors 



It is straightforward to solve the Ward identity (Eq.Q) for the ten longitudinal form 
factors, defined in Eq. ()14|) . in terms of the gluon self-energy function 11 defined by 
n^t(p) = S'^'^ip'^dfiu — PfiPu)^{p'^)- Note that this is not the usual gauge-dependent 
self-energy, but rather the gauge-invariant pinch technique self-energy. At one loop 
in d dimensions this is given by 

(36) 

where Po{d) is given by 

/ ^ 

for massless gluons, quarks, and complex scalars. The mass-dependent results are 
given in section 5 and Appendix E. This result holds for dimensional regularization 
(DREG), whereas for dimensional reduction (DRED) the gluon coefficient changes 
from (7c/- 6) to {8d - 10). 

The longitudinal form factors are given by 





2 


n(p?) 


- ^ipl) 




2 


n(p?) 


- n(pi) 


pI 


-pI 


0, 





(38) 



and of course cyclic permutations yield results for A23, etc. Note that one of the 14 
form factors vanishes to all orders and only the A form factors contain UV divergences. 
In the notation of Eqs. ()20I2H) we have 

AuiQ) = ^j^^liJ^ + J^) (39) 
and similarly for the B and C form factors. 
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3.3 The transverse form factors 



These form factors cannot be determined from the Ward ID, and must be calculated 
explicitly. The algorithm used is briefly described in Appendix A. 

Due to the lengthy expressions, the following shorthand notation will be used for 
the kinematic invariants: 



2 I, 2 



C = pI a=pi-p2 P = P2-P2. l=P3-Pl 

We also define the symmetric invariants 



(40) 



Q = a + /? + 7 
/C = aP + /?7 + 7a 
V = a(3-f. 



(41) 



Note that the dot products can be written in terms of the virtualities a = {c—a — b)/2, 
j3 = {a — b — c)/2, 7 = {b — c — a)/2, or vice versa a = —a — 7, b = —a — [3, c = — /9 — 7, 
but the formulae are simpler and more transparent when selectively written in terms 
of both a, P, 7 and a, b, c. 

We will only write the full gluon contribution explicitly, since the quark and scalar 
parts can be determined from the results of section 3.1 (see Eqs. ()22l23l24p ) and the 
quark-gluon sum rules for the transverse form factors, which are 



;^i2) 



{d-10)( ^ , 2a( Ji - J2) - /3{J2 - J3) - 7(^3 - 



2/C 

(rf- 10) 



/?-7 



J. 



(42) 



The gluon contributions to the transverse form factors in the LT basis are 



+ 



+ 



2K? 



J(lOP + c(/C-7a2-3/37) 



1 - ^^X^) (VJ + a-iJi + a(3J2 + - -^{Ji + J2 + Js] 



7d-6 



2{d~ l){d-2) I 



2p\d-2) 



SV + {d~ 4)ac^ + {d- 2)(4/Ca - c(3-f) 



d-ll 
1 

1 L 



5d-2^ 7d-6 ,^ 
-/C : —ap 



Ji — J2 
a — b 
2d^ - 15d + 14 



d-2 



o 2/j ON 2 7d-6 

2Y{d -2) /C ; 

' ^ ' 2 d-2 



07 



d-2 
2c/2 - I5rf + 14 
d-2 



07 
a[3 



(43) 



{J2 - h) 
{■h - Ji) 
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and 



1 



+ 



2/C2 
d-2 



J 



llC ^{d-2)VQ + {d + l 



abcV 
/C J 



+ V 



d-1 

J + l 



a(/C - 2a7)(Ji - J2) + /3(/C - 2/3a)(J2 - J3) + 7(/C - 2/37) (J3 - Ji 



ci- 1 



- ( Ji + J2 + J3) + 



d-1 



(07 Ji + a/3J2 + jS'yJs) 



(44) 



3.4 The Form Factors in the Physical Basis 

Now we will present the results in the physical ± basis (Eq.©), before symmetriza- 
tion, Eq.p2j). since this is the most convenient way to present the results. Of course, 
these results can be obtained from the relation between the ± basis and the LT basis 
given in Eq.(|16p. but we write them explicitly for future convenience and phenomeno- 
logical applications. 

The quark-gluon sums are given by 



Sqg(^12) 




Sqg(-Bi2) 


4/C V 




{d- 10) 
4/C 


Sqg(-Di2) 


= 




= 




= 0, 



(7 - (3)abJ + {2a + (3)aJi - {2a + 7)6J2 - a(/3 - 7)73) 
■ {acJ + 7 Ji + /9 J2 + CJ3 



(45) 



and the remaining sums (for A23, etc.) are related trivially by permutations a — > 

/9 — ^ 7 — > a, a ^ 6 ^ c ^ a, Ji — > J2 ^ -^3 ^ 
The gluon form factors in d dimensions are 

- 4 /C2 Ai2(G) = afecJ(7/C + /37) + a Ji (7/C/3 + 0^-^ + ^7^^) (46) 
+ h.h (7/C7 + /?72 + ^i^j^ + (7/Ca + P + /Cc-^ ~ ^ ^ 



d-\ 



.2a(7rf- 6) + 7 



-4/C2 5i2(G') =a6J(7/C + /37)(7-/3)+aJi(7/C/5-67(/5-7) + /C 

- 6J2(7X:7 + a/3(/3 - 7) + ^ M7^-6)+/^ j ^ _ ^) j3(^7^^ + p + ^^"^ " ^ 



d-1 
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16 Ci2{G) = cJ 3/C^(10a + c) + /C(a^ - 60/^7) + VlC{d + 4) - V{d + l){a^ + 2/37) 



+ iK 



(3 - 2rf) _ 30d + 24) 



d- 1 



-a — 7- 



d-l 



3/3 



- — -+7(2Q-3a) 



d-l 



d-l 



d-l 



P'y - 47 



(3 -2d) (d'^-30d + 24) 

■a — p — h 37 



d-l 



d-l 



2 (^-2) 
d-l 

V{d+1) 



(47) 



c/- 1 



d-l 



( 30rf-31) ^, ^ ^ 
d — 1 



- 4c 



rf- 1 ' ' ' d-l 
{d-2) (d^-Ad+l) 



d-l 



d-l 



/37 + l)P(2Q-3a) 



16 /C^ Di2(G') = a6(a - 6) j(/C(Q + 2a) -{d+ l)p) - aJi 



(48) 



/C 



d-l 



d-l 



d-l 



{2a + p)V{d+l) 



bJo 



id' - 4) ^2 



(rf2-3) (c/2-4c/+l) 2(^^-2) 



d-l 

- {2a + -f)V{d+l)\ +{a-b)Js 



d-l 

(4rf - 7) 
d-l 



— 07 



d-l 



4a' 



d-l 



3a2 - /37 



(d^ - 3) 
d-l 



-aP(d+ 1) 



16IC^H{G) =abcj(v{d+l) - ICQ 
+ aJi 



+ 6J2 
+ CJ3 



3- 


- 2d 


d 


- 1 


3 - 


-2d 


d 


- 1 


3 - 


- 2d 


d 


- 1 



d2 


-3 


d 


- 1 


d2 


-3 


d 




- 1 


d2 


-3 
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■a7-/32j/C + (d+l)/3P 
a/3-72 /C + (d+l)7P 
/37-a2l/C + (d + l)aP 



(49) 



16}C^S{G) = {a-b){b-c){c-a)j(3}CQ-{d+l)V 



+ (&-c)Ji| |^ + /C 
+ (c-a)J2| |^ + /C 
+ {a-b)J,\^ + }C 



Ac' 



,d-2 
d-1 

d-2 
d- 1 

d-2 
d- 1 



07 



d^ 


-4d4 


- 1 




rf- 1 




d^ 




^ 1 




d-\ 








1 


d-\ 



(50) 

- (d+l)P(2Q-3/3) 
-(rf+l)P(2Q-37) 
-(d+l)P(2Q-3«) 



Now we turn to the physical symmetrized basis. From Eq. ()12|l . we see that for 
any triplet of form factors, say A^^ we have 



A_ 

where we have defined 



- ( A12 + A23 + A31 



Ax2 + AA23 + 

AA31 



- ( A12 + AA23 



+ ZA2 



Ax - tA2, 



(51) 



^1 
A2 



^(^12-^(^3 + ^1; 
6 V 



(52) 



and A2 correspond to the real and imaginary parts of A± only when J, Ji, J2, J3 
are real. This occurs (in the massless case) when /C > 0, which can only happen if all 
three gluon virtualities are of the same sign, either all spacelike or all timelike. This 
is often not the case for real problems. In general, however, it can be shown that 



A*^{a, b, c) 
B*{a,b, c) 
Bl{a,b,c) 
Coia,b, c) 
C![.(a, b, c) 
Do{a,b, c) 
D^{a, b, c) 
H*{a,b, c) 
S*(a, b, c) 



A:f{—a, —b, — c) 
So(-a, -6, -c) 
i?=p(— a, —b, — c) 
-Co{-a, -6, -c) 
— C=p(— a, —6, — c) 
-i:)o(-a, -b, -c) 
— Z^=p(— a, —6, — c) 
—H{—a, —b, — c) 
— 5'(— a, —6, — c). 



(53) 
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Furthermore, all of the above form factors except for Aq are scale invariant, 

F{\a,\b,\c) = F{a,b,c) = F{a/c,b/c,l) A > 0. (54) 

Only Aq is not scale invariant and does not satisfy a simple reality condition. 
The quark-gluon sums are given by 

SQG(i?±) = ^^^^^|^-3(Q$2 + /CJ(/?-7))±zv^(Q$i-/CJ(Q-3a^ 

Sqg(C±) = ^-^^^(^^1^'^^^^ (55) 

SqgPo) = SqgP±) = 
Eqg{H) = Sqg(5)=0 

where we have defined the commonly occuring functions 

$0 = CLbcJ + apji + 67J2 + caJs 

$1 = (/C-3/37)J-37Ji -3/3J2 + 3(/3 + 7)J3 (56) 
$2 = a(/5-7)J+ (2a + 7)Ji - (2a + /3)J2 + (/3-7)^3. 

From the definition of Sqg in Eq. ()22|) we see that the quark and gluon (and thus 
scalar, by Eq. ()23p ) contributions have the same functional form for the seven form 
factors which have a zero in the above. Letting F stand for A±, Dq, D±, H or S, we 
find that 

F = {n^-^Y. ^f^f + dh^ ^'^') ^^^^ ' ^^^^ 

/ 

which, in d = 4 QCD, reduces to F = ig'^{N^ - Nf)F{G). 

In addition, both Aq and Bq are governed by one function, since they satisfy 
different sum rules. In particular, the tree level tensor structure has coefficient 

A,^-%(11c.-^J^Tt,N,-^^J:t.N.]%. (58) 
" 21c y 3 3{d-2) ^ ^ ^ 3(ci-2)^ ' y " ^ ' 

This form factor will be discussed in more detail in the next section. 
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Also, one finds from explicit calculation that the scalar contribution to Bq van- 
ishes, 

BoiS) = BoiG) + 45o(g) = 0, (59) 



and thus 



Bo = tg^ 



[nc-\Y, TfN^ Bo{G) where 



Bo{G) = -^[{a-b){b-c){c-a)J+{b-c){2Q-3f3)Ji + {c-a){2Q-3j)J2 
+ {a-b){2Q-3a)J3j (60) 

Finally, since S = exactly, we know that our fourteen dimensional basis is 
degenerate, which is reflected in the fact that S + SDq = 0. Hence we define a new 
basis tensor (Iq = do — 3s so that -00*^0 + Ss = Do{do — 3s) = Dod'^. 

Thus, we find that eight of the thirteen nonzero form factors have the same func- 
tional form for gluons, quarks, and scalars. Only the five form factors B±, Cq and C± 
do not. These statements are basis dependent. One can always find bases where none 
of the form factors have a vanishing QG sum rule. In the course of our calculations, 
we found that the ± basis gives the maximum number of such zeroes among bases 
which are reasonable and contain the tree-level tensor structure. In this sense the 
(symmetric) ± basis is the simplest and most compelling. We will see in the next 
section that this is also the most convenient basis for perturbative calculations. Of 
course, as discussed in a previous section, with supersymmetry every form factor is 
proportional to j3o, and so supermultiplets are governed by the same function in any 
basis. 



4 Three-Gluon Vertex in Perturbation Theory 

Applying the pinch-technique (PT) construction to the three-gluon vertex occurring 
in a physical process involving three external on-shell legs, one arrives at a dressed 
tree-level skeleton graph, dressed with pinch-technique vertices and self-energies as 
shown in Fig. 2. Generically the amplitude of the three-gluon graph can be written 

M = C(7oVi'^i\/2'^^\/3'^^D^,,,(A;i)D^,,,(A;2)I^;.3.3(fc3)r^,M2M3(^i, h, ks), (61) 

where C is the overall color factor and go is the bare coupling. The PT vertices Vi are 
for gluons coupled to external particles, whose indices are suppressed; this is shown in 
Fig. 2 for external quarks. F = Yli=i ^ifi gauge invariant three gluon vertex. 
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(A) 



(B) 
+ pemis 



Figure 2: The tree-level skeleton graphs dressed with pinch-technique vertices and 
self-energies, which are used to define effective charges both for the quark-quark-gluon 
vertex (depending on a single gluon virtuality) and for the three-gluon vertex, which 
depends on three different momenta. 



whose thirteen form factors (Fj) are given in the preceding section. Finally the "gauge 
invariant" PT gluon propagator is 



k,,h 



, / r \ Quv TT- ■ covariant ■ 

VW = |jLn4^|-| axial ^g^^g^^ 



l^k) = I I in I I gauges, (62) 



where ^ is the gauge fixing parameter. This is "gauge invariant" in the maximal sense, 
i.e. the gauge dependence comes only from the tree level terms, and in particular 
n(A;2) is totally gauge invariant. 

Regardless of whether the external particles are quarks, gluons, or scalars, the 
vertices satisfy Vi,^ = when these particles are on shell (OS). One can then show 
that the gauge dependent terms coming from Eq. ()62|1 vanish in the full amplitude 
consisting of all of the graphs in Fig.Q. This can be seen trivially in the covariant 
gauges where the gauge cancelations occur graph by graph, and with some work in 
axial gauges, where the cancelation occurs between all of the graphs. In the latter 
case, one must use the fact that the three gluon vertex satisfies the Ward ID in Eq.©. 
Therefore we can take 

^'■"^^^ ^ F(i + n(P))- 
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Also, in the ± basis (Eq.((7j)) any tensor with a in any slot gives vanishing con- 
tribution to ^A. For example, (+00) = + k^j^j^g^^f^a = ~kl,^^l9^^2^J.3^ ^i.^ti dots 

into Vi^^i yielding zero. Hence only (00—), (—00), (0—0), and ( ) contribute, and 

we find 

^ = {i+uikim+nUm+uiki)) \ikiki fcs), (64) 

where the three-gluon vertex connected to on-shell (OS) external particles is 

r?fM2M3(^i' ^3) = (1 + ^)«o + A+a+ + A„a„ + Hk (65) 
in the notation of section 2 where hatted objects are three index basis tensors. 
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Figure 3: The effective scale Qgjj(10 GeV^, 10 GeV^,p^) is the lower blue curve, 
while Qe//(-10 GeV^-10 GeV^p2) = g^^/lO GeV^ 10 GeV^ -p^) is the upper 
red curve. These both asymptote to zero, although very slowly for the upper curve. 

Now one naturally defines PT effective charges by^ 



Since we only have a single power of Qq for each l/(l+n(A;j )), this leaves Y\i=i V V 1 + n(/i;, 
1 — |(n(/i;^) + n(/c2) + n(/c|)) to be absorbed into the three gluon vertex. Thus we 
have 

3 1 

n^^^7TfTr™(^i'^2,/c3) = {l+A,)ao + A+a+ + A_d_ + Hh 

Mkl kl kl) ^ Aoikl kl kl) - ]^{li{kl) + Yi{kl) + Ii{kl)) (67) 



•^Ea. H66|) holds for external fermions or scalars. but for gluons one would instead have three 
additional three gluon effective couplings, as is clear from the derivation of Ea. (|69|l . 
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1 2 4 6810 20 50 100 200 400 
p'*2(GeV'^2) 



Figure 4: The effective scale (5e//(10 GeV^,p^,p^) is the lower blue curve, while 
(5e//(-10 GeV^,p^,p^) = (5e//(10 GeV^, -p^, -p^) is the upper red curve. These 
both asymptote to lOe^ GeV^ 31.25 GeV^. 



and 



"^1 "^2 "'S 



-^0 



;i + Ao)do + A+a+ + A„a_ + Hh . (6. 



This naturally leads to the effective coupling of the three gluon vertex 



g{ki, ^3) 
5(a, b, c) 



goil + Aoiklklkl)) 
]p{a,b,c) ao 



47r l-2Ao(a,6,c) 



(69) 



first obtained by Lu in jSH]- Our amplitude then takes the final form 



M=Cg{kl)gikl)g{klMklklkl 



ao + A+a+ + A_a_ + Hh . (70) 



Recall from the previous section that A±, H cc Nc — Nf in QCD. 
The three-gluon effective coupling evolves according to 



a{a, b, c) 



a(ao,6o, Co) 



1 -2[Ao{a,b, c) - v4o(ao,&o,Co) 



(71) 
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In four dimensions with regularization scheme R = DRED or DREG we have 

Ao{a,b,c) = --^po L{a,b,c) -log - Cuv - Vs 
on I 

11 2 1 
where /5o = ^N, - -Nf - -N, 

6 6 

Cuv = ^-7i? + log47r (72) 

r]3 = {2 + n) + —^Sr^dreg 

Q = ^Cl2(7r/3) ^ 3.125 

The scheme dependence Sr^dreg is explained in more detail in Appendix D. Here we 
have defined 

L(a, ^^^) ~ ~jQ {^'^ + b + log c — abcJ{a, b, c) j + Q, (73) 

and the (massless) triangle integral function J = J{a, b, c) = — IGzvr^ J(a, b, c) is given 
in Appendix B in terms of Clausen functions, Eqs. (|121ll22p . This result fEqs. (l69|72|) ) 
differs from Lu jSH] by only the finite constants which (slightly) affects the numerical 
extraction from data. The discrepancy can be traced to the inconsistent application 
of dimensional regularization in [38] . 

The logarithm-like function L satisfies 

L(a, a, a) = logo (74) 

since 

J(a, a, a) = ^a2(7r/3). (75) 
av3 

One can use the real part of this function to define an effective scale of the three-gluon 
vertex: 

L(a,6, c) = log ((5e//(a, 6, c)) + ilm L(a, 6, c) 
Qlff{a,b,c) = |ar^/'^|6r^/'^|c|^^/'=exp(^fi-^Re7(a,6,c)). (76) 

This is sensible since the dimensions of Qljj-{a, b, c) are indeed mass squared. 

The three gluon effective charge a{a, b, c) is related to the usual MS coupling 
a{q^) by 

1 ... ^2 



a{a,b,c) a{iJ?) Atx\ jji 



a{e-'''Qlff{a,b,c)) 47r 
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+ i^lm L{a,b, c). (77) 




Figure 5: A contour plot of Qljj-{l,x,y). The contours, from red to blue, are at 
2,4,6,8,10,12,14,16,18,20. 



Since exp (^) ~ 14, we see that when using MS, the scale should be fourteen times 
lower than the typical virtualities of the gluons, given by Qeff{a,b,c). Of course, 
this is true only if the three-gluon vertex diagram dominates the physical process. In 
general there will be different scales at the various quark-gluon vertices when using 
the PT scheme (as seen in Eg. (170])). In contrast, in MS the same scale is used 
at every vertex. The following approximate values of the three-gluon coupling are 
derived from Eq. (|77jl (including the effects of quark masses which discussed in the 
next section) for various symmetric timelike(T) and spacelike(S) configurations : 



SSS : 5(-M|, -M|, -M|) 0.192 

f|, -Ml +M| 



SST: 5(-M^,-M|,+M^) « 0.157 + 0.023/ 



STT: 5(-M|,+M|,+M^) ^ 0.156 + 0.025/ 



r2 I /\#2 I /\/r2^ 



TTT: 5(+Mi,+M|,+M|) ^ 0.170 + 0.062/ 



It is clear that the three-gluon coupling is stronger than naively expected from 
aMsiMz)^ 0.118. 
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Figure 6: A three-dimensional contour plot of Qlfjr{l,x,y). 
The effective scale QlfAa,h,c) satisfies the following relations: 



Qlff{a,b, c) 


= Qlffi- 


a, —6, — 


c) 


Qlff{Xa, Xb, Ac) 


= IMQlff 


(a, 6, c) 




Qlffia, a, a) 


= \a\ 






Qlff{a, -a, -a) 


^ 5.54 a 






Qlffia, a, c) 


^ |c|e^-2 


for 


\a\ ^ c 


Qlffia, -a,c) 


^ c 


for a 


> kl 


Qlffia,b,c) 


a 


for 


|a| > |c| 



Lu [SHI previously found the last of these limits in the case where all momenta are 
spacelike, giving an effective scale QminQmed/ Qmax- It should be noted that the rate of 
convergence to the above limits strongly depends on the signatures (S = spacelike 

< 0, T = timelike <-> > 0) of the virtualities a, 6, c. If the signatures are mixed 
(TTS) or (TSS) then the convergence is very slow, and the effective scale tends to 
stay larger compared to the cases (SSS) or (TTT). 

Some plots demonstrating the novel behavior of Qljj are given in Figs. ()3l4l5l6|) . 

5 Phenomenological Effects of Internal Masses 

So far, all fields propagating in the triangle graphs have been treated as massless. This 
was useful for simplifying the discussion and elucidating the general structure of the 
radiative corrections and the = 4 sum rules. However, in real world applications 
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one usually does not have all three gluon virtualities in the same desert region Mj -C 
a,b,c <^ Mj+i. Thus, mass corrections should be taken into account. We have 
calculated the effects of massive fermions (MQ), massive scalars (MS), and massive 
gauge bosons (MG) for all of the form factors; the complete results are given in 
Appendix E. The corrections for the case of massive fermions were first obtained in 
Ref. |2ni and we are in agreement. Here we will focus on the massive quark (MQ) 
contribution to the form factor multiplying the tree level tensor structure, which from 
Appendix E and section 3 is 

^j^2 2^ 8 r 

Ao{MQ) = _ ^^ Jm + 5^(^_ 2) ^^'^'^M + apJiM + &7<^2Af + caJ^M ■ (79) 

Here Jm, Jim, J2M, and J-sm are the massive analogs of J, Ji, J2 and J3, respectively. 
The two-point function Jim and tadpole Tm are reviewed in Appendix A, while 
Appendix B is devoted to a discussion of the massive triangle integral, Jm, and its 
analytic continuations and various limits. 




2 4 6 8 10 12 14 16 
a/M''2 



Figure 7: LMqia/M'^, a/M^, a/W^) vs. a/M^ for timelike a > 0. The sohd line is the 
real part and the dashed line is the imaginary part. 

As in the previous section, when considering a physical matrix element we always 
have the combination Aq = Aq — ^ (Hi + H2 + H3) multiplying the tree-level tensor 

structure. This leads us to consider the massive quark contribution Aq{MQ) = 
Aq{MQ) - |(Hi(Mg) +H2(Mg) +H3(MQ)), which upon using Ea. (inH|) . inserting 
the prefactor ig"^, and expanding around d = 4 becomes 



Ao{MQ) 



+ 



<y.s 
2M2 



uv 



log 



1^ 



2 
3 



3/C 



[abcJM — a'yC{a) — a(3C{h) — /3'yC{ 



J, 



AI 



2 - C{a) 2 - C{b) 2 - C{c) 
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Figure 8: Lmq{-QVM^, -Q^M^ -Q^M^) vs. Q/M for spacelike -Q^ < 0. 



Here J, 



M 



— IGin'^JM and 



C{a) = f (a) loj 



via) 



via) - 1 



where v(a) 



4(M2 - ze) 



1 



comes from Jim and has the analytic continuations given in Eq. p06|) . The three-scale 
logarithm-like function for massive quarks (MQ) is thus given by 



fa b c \ 



AoiMQ) 



— ( a7£(a) + af3C{b) + f3'yC{c) — abcJM{a, b,c)) +Q 



2M^ 



C(a) 



^c) - 2 _ 
'Jm 



(82) 



Vl27r/ 



This massive logarithm-like function has the following limits : 

a b c 



T ( ^ A- ^\ 



2 + n > |a|, |c| 

L(a,6,c) - logM^ < |a|, |c| 



^3) 



with the number Q f» 3.125 defined in Eq. ()72|) . The convergence to the massless limit 
is very slow, indicating that threshold effects must be included for most applications. 

In Figs. ()7p8|l we have plotted Lmq for the symmetric case a = 5 = c for time- 
like and spacelike momenta. For the timelike case, the threshold at a = 4M'^ and 
the pseudo-threshold at a = 3M2 are evident. In Fig.Q the mixed case (TTS) 
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a/M"2 



Figure 9: LMQ{a/M^,a/M^,-a/M^) vs. a/M^. The solid line is the real part and 
the dashed line is the imaginary part. The real thresholds are at a = ±4M^ while 
the pseudo-threshold is at a = 5M^. 



LMQ{a/M'^,a/M'^,—a/M'^) is plotted, where the thresholds at a = ±4M^ and the 
pseudo-threshold at a = 5M^ are evident. For the purely timelike (TTT) case in 
Fig.((Zj), there is a discontinuity in the imaginary part and the real part diverges at 
the pseudo-threshold. In contrast, for the mixed signature case (TTS) the imaginary 
part diverges while the real part is discontinuous. This pseudo-threshold phenomena 
is explained in more detail in Appendix B. 

From the above results, one can define the effective number of active quarks which 
characterizes the effects of quark mass : 



a b c \ 
M2' M2' M2y 



^ L 



a b c 



B4) 



This clearly goes to zero and one in the limits JVP ^ |a|, |c| and <C |a|, |c|, 
respectively. 

To motivate this definition, let's look at the single-scale pinch-technique effective 
charge (using the notation of as a function of spacelike momenta a = —Q^ 



(85) 



where (3p is the contribution of each particle p to the first (3 function coefficient, and 
to a very good numerical approximation (for spacelike momenta) 



LpiQ'/ml) ^ log 



mi 



32 



10 
Q/m 



Figure 10: The effective number of quark flavors. The lower blue solid curve 
is iVi?(Q^/M^) for the symmetric spacelike {a = b = c = —Q^) three-gluon vertex, 
while the upper dashed red curve is the fermion number of flavors Ni/2{Q'^ /M"^) for 
the single scale effective charge. 



where the constants r]p are 5/3,8/3, and 40/21 for massive fermions, scalars, and 
gauge bosons, respectively. The exact one- loop formula are given in Eq. (23-26) of 
Ref. |13j. although the analytic continuation in Eq.(26) of that paper should have 
opposite imaginary part. This effective charge satisfies the RGE 

da{Q') ^ dL,{Q^/ml) _ ^ f 

(ilogQ2 ^71^^" dlogQ^ ~ All ^^" "^ [ml 



g2\ _ dL^iQ^ml) dLpiQ^m. 



2^ 



"ymlJ dlogQ^ dhgml ^ i + ^e^.' ^ ^ 

The function Np goes to one when Q"^ ^ and zero when ^ and unambigu- 
ously measures what fraction of particle p is "turned on" at scale Q^. 

Moving back to the three-scale case, we now have the complication that our ef- 
fective charge is a solution of a multi-scale RGE 

d5(a, 6, c) _ sr^ n d f _a b 

dloga ~ 47r^^^dloga '''^ V^^'M^'M^J ' ^^^^ 

and two other permutations with a — s> 6 or a — >■ c. This leads to three different Nf's 



N 



I' a b c \ d fa b c \ 

M2 ' M2 J " rfloga M2 ' 



^9) 
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and two cyclic permutations, each of which goes to 1/3 in the symmetric desert 
a = b = c ^ M^. This suggests adding all three together to define a symmetric 



at(^^^\('^ d d \ I a b c \ 

which is in fact the same as given in Eq.(jHH). 

The results for Np can be obtained with the help of the results in Appendix 
B. Instead of presenting these lengthy results, let us focus on the symmetric case 
a = b = c, where we find for spacelike a < 

/a\_ d fa a a \ d fa a a \ 

^ \Jp) " dl^ ^'^^ ' M2 ' Jp) ^ ~d\ogM^ Im2 ' M2' M2 J 

= 1 + 18- + 2MV„(a, a, a) + ,iM^C(a)^^-^^^^,^^^^-^91) 

In this example, spacelike a is chosen to avoid the pseudo-threshold a = 3M^ and 
the threshold a = 4M^. Fig. ljlUj) shows a plot of this along with single-scale quark 
number of fiavors function Ns=i/2 from Eq. (j87|) . 

The negative value of Np at ^ Q ^ AM is not entirely novel, as a similar 
behavior was also found in the context of two-loop quark mass corrections to V- 
scheme effective charge jlU]. It is essentially due to the anti-screening of color charge, 
this case in the triangle interaction, and does not arise in the one-loop single scale 
effective charge, as seen in Fig.Q. 

Using the results of Appendix E, the above analysis can be easily extended to the 
case of massive scalars (MS) or massive gauge bosons (MG), which have the following 
logarithm-like functions 

^^^^ (i^' ii^' " ^(a-fC{a) + apC{b) + f3-fC{c)-abcJM{a,b,c)^+n 

- + ^4 + Jj,i (92) 



^^^{^'JP'^J = ^(a-fC{a) + aPC{b) + p-fC{c)-abcJMia,b,c)^+n 



+ ^,.fM^ + ^+A£)^_7;^|, (93) 



7 

The qualitative behavior is the same as the quark case. 

Finally, we should consider the limitations of the effective scale Qljrj{a,b,c) in- 
troduced in the last section and effective number of fiavors NF{a/M'^,h/M'^,c/M^) 
discussed in this section. Given the complicated structure of the full mass dependent 
form factors, such tools for characterizing and understanding the behavior of the ver- 
tex are helpful. However, in a real calculation such methods may be of limited use 
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and the full mass dependent results should be used. For example, the effective scale 
Qeff has been defined only in the massless case so far because the definition becomes 
complicated and somewhat arbitrary in the massive case. In particular, consider the 
possible definition (for QCD): 



where A''^ is some suitably defined number of flavors, possibly a step function such 
as 6{a + b + c — 3Mg), possibly the Np defined in Eq. ()84j) . or some other definition. 

It should be clear that any choice of Nq determines Qlff, and vice versa, and there 
seems to be no compelling choice for these quantities. Furthermore, in the approach 
advocated here, the couplings at each vertex depend on physical momentum scales 
which will typically be integrated over in the phase space. Thus, matching onto a 
conventional MS type approach can only be done at the end of the calculation, so 
that trying to define a Qljf at an intermediate stage is not very useful. 

Thus, in real world applications, one should generally use the full results for the 
mass dependent form factors. This constitutes a multi-scale analytic renormalization 
scheme that contains information which cannot be obtained in the simple single- 
scale leading-log renormalization methods. In other words, every three-gluon vertex 
(at tree-level) can be dressed, or "RG improved" , with this gauge-invariant effective 
coupling and the associated form factors, which are process independent and contain 
more information than the MS procedure. 

6 Conclusions and Future Directions 

The results of this paper represent only a fraction of what is needed for a re-organization 
of perturbation theory into fully gauge-invariant pieces with physical content, each of 
which can be renormalized independently, leading naturally to a physical multi-scale 
analytic renormalization scheme. This is possible due to the remarkable properties 
of the pinch technique (PT)/ Background Field Method quantum Feynman gauge 
(BFMFG) Green's functions. There is still much progress that can be made in cal- 
culating these Green's functions in perturbation theory. 

The present paper gives a complete and general characterization of the off-shell 
three-gluon vertex at one- loop. A similar study of the gauge-invariant triple gauge 
boson vertices of the Standard Model jHj would be very useful. It may also be possible 
to quantitatively look at the unification of triple gauge boson vertices and couplings, 
in analogy with the work on the unification of single-scale PT couplings^^l- Some 




+ log M 

Ml' M^' Ml) ^'"S'""" 



) 



(94) 



35 



progress has been made on the conventional gauge-dependent three-gluon vertex at 
two loops 1151121; which gives hope for eventually treating the gauge- invariant three- 
gluon vertex at two-loops. 

The gauge boson two-point functions and the associated effective charges for QCD 
electroweak theory ^HIIHIC^; and supersymmetric grand unified models ^3] 
have been calculated in the past to one-loop. We know from general principles the 
divergent parts at two-loops, but no complete two-loop PT calculation as yet exists. 
To fill this gap, the two-loop gluon PT self-energy will be presented in the near future 
[HZj. This will allow for a more precise determination of coupling from data, as well 
as giving the three-loop (3 function coefficient. Furthermore, by the Ward identity in 
Eq.Q, this also yields the longitudinal form factors of P^^i^aVa through two loops. 

The gauge- invariant PT/BFMFG quark self-energy turns out to be equal to the 
conventional self-energy in the Feynman gauge [2^], and so is known through two- 
loops US]. Due to the Ward identity |211 satisfied by the PT/BFMFG quark-gluon 
vertex, this also yields the longitudinal form factors of F^^'^^^ through two loops. 

In QCD, another logical step is the four-gluon vertex at one-loop. In the general 
off-shell case, there are hundreds of independent tensors structures and form factors. 

Beyond perturbation theory, the study of Dyson-Schwinger Equations [3] and 
renormalons [Zj in the PT/BFMFG approach may yield new insight. 

To summarize, in this paper we have analyzed the behavior of the gauge-invariant 
three-gluon vertex at one-loop. Starting from the symmetry principles governing 
the vertex, a convenient tensor basis decomposition was given in Eqs. fl7ll2|) . As 
seen in Eq. ()65|) and the subsequent discussion, this basis is the most convenient for 
phenomenological studies, since it is built out of "transverse" (— ) and "longitudinal" 
(-I-) momenta, the latter of which vanish when dotted into external on-shell vertices, 
thus leading to relatively simple matrix elements. In the case considered in section 
4, only four form factors remain, rather than the thirteen which would be present 
in a generic basis. Nonetheless, the choice of basis is only a matter of convenience, 
and the real physics lies in the thirteen non- vanishing form factors given explicitly in 
section 3. 

The supersymmetric relations between the scalar, quark, and gluon contributions 
leads to a simple presentation of the results for a generic (unbroken) gauge theory. 
Only the gluon contributions to the form factors are given explicitly in section 3, 
while the quark and scalar contributions are inferred from the homogeneous relation 
Fg + 4Fq + (10 — (i)F5 = and the results for the relatively simple sums Sqg(-F) = 
(d-2) _|_ "vvhich are given in section 3 for each form factor F. The extension to 
the case of internal masses is outlined in Appendix E and leads to the modified sum 
rule Fmg + 4Fmq + (9 - d)FMs = 0. 

The phenomenology is largely determined by the form factor of the tree-level 
tensor structure, which in section 4 is used to define a three-scale effective charge 
5(a,6, c). In addition, the characteristic scale QlfAa.h.c) governing the behavior of 
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the vertex and the effective charge was analyzed, thus providing a natural extension of 
BLM scale setting to the three-gluon vertex. Physical momentum scales always 
set the scale of the coupling. The phenomenological effects of quark masses are 
discussed in section 5 and are found to be important for generic physical applications, 
since decoupling is slow and a complicated threshold and pseudo-threshold behavior 
is observed. An important next step is to fully apply these techniques to a physical 
process. In the future we will present such an analysis for the hadronic production of 
heavy quarks, where the importance of the form factors other than the tree-level one 
(Aq) will be addressed. The interpretation of the pseudo-threshold phenomena also 
deserves further study. 

Acknowledgements M.B. would like to thank Lance Dixon for useful dis- 
cussions regarding the second order formalism of the BFM and the supersymmetric 
decomposition of one-loop amplitudes. 
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Appendix A: Reduction to Scalar Integrals 



First, we will describe the evaluation of the massless integrals, and then briefly 
mention the modifications due to internal masses. As before, we will use the shorthand 
notation 

a = pl b = pl c = pI a = pi-p2 P = P2-P3 1 = P3-pi (95) 

In order to evaluate the integrals in an efficient manner, it is very convenient to 
choose a manifestly symmetric routing of the loop momenta, as shown in Fig^ where 
clearly 

h=P2 + h k=P3 + h h=Pi + h- (96) 

Of course there is only one integration momenta /, which can be chosen to be /i, or 
/s, thus breaking the cyclic symmetry. However, using the symmetric labeling greatly 
simplifies the analysis. 

First we decompose the full vertex F into longitudinal (L) and transverse (T) parts, 
T = Tl + Tti as in Eq. ()14|l . The tensor integrals in Eq. ()18|l are then converted into 
scalar integrals by applying projection operators. In doing so, the longitudinal (L) 
and transverse (T) parts essentially decouple, and the ten independent L form factors 
are easily found either directly, or by solving the Ward ID, resulting in Eq.()38p. The 
remaining four T parts are found by applying the following four projection operators 
to Eq.dlHI) : 200,030,001, and 231, where as in Table 1 we have defined 030 = 
P3^J.29^J.l^J.3^ etc. Thus, for the gluon contribution G we have four scalar integrals : 
G(200) = p2/xi5'/x2M3^MiM2M3' ^(030), (7(001), and (7(231). Similarly, there are four 
integrals for the quarks and scalars as well. In the numerator of each of these integrals 
there will be various dot products of momenta, which can always be reduced to 
momentum squares using Eq. (|nH|l . For example, pi-h = (/| — /| —p\)/2 and h-l^ = 
(^1 + ^3 ~ Pi)/^- Thus we are left with integrals of the form 



/2/2/2 ' V^'/ 



where J = J -(0^ and i,j,k E {0,1,2}. Using the standard rules of dimensional 
regularization, it is easy to see that any integral with any two of i,j, k nonzero must 
vanish. Furthermore, it is straightforward to show that 



'200 



—Phoo I020 — —ihio ho2 — — tt-^ooi- (98) 



Thus we are left with only two types of integrals: (1) the trivial two point integrals 
Ji, J2, and J3, where 



- ~ ' wi'j ni+PiY' ^^^^ 
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and (2) the master triangle integral 



J = J{pI,pI,pI) = f 



(100) 



For the gluon contribution, for example, one then has a system of four equations 
with four unknowns, the transverse form factors. Denoting the gluon contribution to 
the longitudinal projections by Lg'(200) = 200Ti(G), etc. we solve for the transverse 
form factors 



FjiiG) 
\ H{G) ) 



M- 



where M' 



T 



-/C 



/ ^(200) 


- ^g(200) \ 










G(030) 


- ^g(030) 










G(OOl) 


-^g(OOI) 










\ G(231) 


-^g(231) ) 










( ^ 


{d-l)p 




/9 


2 


-d\ 


7 


7 {d 


-1)7 


2 


-d 


{d-l)a 


a 




a 


2 


-d 




/C 




/C 




) 



(101) 



and similarly for the quark and scalar contributions. 

The above procedure can also be followed for the massive case, with only a few 
modifications. First, the tadpole M'^Tm = J jr^j^ does not vanish. Thus, instead of 

J ]Tp = ~PJi we now have 



{II 



M2) 



M, 



where J 



IM 



(P-M2)((i+pi)2-M2)- 



(102) 

We also need the following 



result and permutations : 



q - 



q - 

q-M^ 



aM'^T, 



M- 



Finally, we have the master triangle integral with nonzero masses 

1 



J, 



M 



Jm{Pi,P2:Pz) 



{q- M^){q- M^){q~ M^) 



(103) 

(104) 
In the 



To summarize, in the massive case we need Jm, Jim, J2M, J3M, and Tm- 
massless case we need J, Ji, J2, and J3. For each of these we pull out the factor 
and define Jm = j^Jm, etc. 

Some formula for these integrals in d dimensions can be found in jHS]- Here we 
will give only the expansions in four dimensions and define Cuv = j ~ Ie + log47r 
where ci = 4 — 2e. 
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The tadpole integral is 



Tm = Cuv + I -\og—. (105) 



The two point integral is 

Jim = Cuv + 2 -log— -C{a) (106) 

r 2i;tanh"^ (t;-i) = t;log^ ^ f a<0 

£(a) = J 2Utan^^ (U-^) " I for < < a < 

[ 2i; tanh"-^ (v) - iirv = v log - zvrt; J [ a> AM^ 

and the generalized velocities are 



; 4M2 /4M2 , 

y = Wl u=W 1. 107 

a \ a 



In the massless limit this becomes 



Ji = Cuv + 2- log ^ + iTTe{a). (108) 



Appendix B: Results for the Triangle Integral 

The massive triangle integral 

/d^l 1 
{2nr {If - + ^e){q - + ze)(/i - + ze) ^^^^^ 

is finite in four dimensions. We will give the results for Jm = — zlGvr^ Jj\/- This integral 
has been discussed previously in the literature [IS] jlHI |1ZI (Ml • In particular, 'tHooft 
and Veltman [33] derived a formula which is valid for all values of the kinematic vari- 
ables a, b, c and mass M, although careful analytic continuation is required. We will 
first write the results of [IHj in our notation and then discuss the analytic continu- 
ations. The various functions involved and some reference formula are summarized 
below in Appendix C. 

Defining p = y/—IC, where as before 

fC = aP + + -fa = ~ (a^ + + - 2{ab + bc + ca)) , (110) 
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we have 

J^i = -^iHa\b, c) + hib\c, a) + I-sicla, b)). (Ill) 
zp 



The results for l3{a\b, c) can be expressed in terms of the velocity v = y 1 — ^^^"^^"^"^^ 
and the variable x = f3/p : 

'Z2 

l3{a\b,c) = Li2{zi) -Li2{zi) + Li2{z2) -Li2{z2) - r]{x - v,x + v)log — 

+ ri{-l-v, ) log 2:1 -77(1 -w, )\ogz^ 

X — V X ~ V 

+ 'q{-l + v,—^)\ogZ2-'n{l + v,—^)\ogzl (112) 

X + V X + V 

where we have defined 

x + l x + l _ x-l _ x-l 

Zl = Z2 = ■ Zi = ■ Z2 = (113) 

x — v x + v x + v x — v 

and the function ri{x,y) compensates for the branch cut in the logarithms: 

logxy = log X + \ogy + rj{x,y) (114) 

rj{x^y) = 2?7r^0(— Im x)^(— Im y)6'(Im xlm y) — ^(Im x)^(Im ?/)6'(— Im xlm y) 

The other two integrals l3{b\c, a) and /3(c|a, b) are easily obtained by permutation 

of the above results, so that x = 'y / p,v = y^l — and x = a/p,v = \j ~ re- 
spectively. Although these results entirely characterize the massive triangle function, 
it is a rather tedious exercise to analytically continue the results to the six different 
physical kinematical regions. To our knowledge, such complete analytic continuations 
have not appeared in the literature thus far. 

Jm takes different forms for /C > and /C < since then the variable x is 
imaginary and real, respectively. The case /C > can occur only if all momenta 
are spacelike (a, 6, c < 0) or timelike (a, 6, c > 0). The case )C < can occur for 
momenta of any signature. Thus, if all momenta are spacelike or all timelike, the 
ratios of momenta will determine if /C > or /C < 0. For each of these two cases, 
we must also distinguish when a is spacelike, timelike below threshold, and timelike 
above threshold. For timelike above threshold and spacelike, the generalized velocity 

1 — is real, except for the ie term which is used in the analytic continuation 



and hence not included below. Below threshold v = iv = iy — 1. 

Case JC > 

For /C > we have 

p = V— /C = iVK. = ip 

X = (3/p = —1(3/]) = —iw (115) 
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/C > and v real <^=^ (a < or a > AM'^) 



iM^.c) = d 2Cl2(20i)-Cl2(20i- 202) -012(201 + 202) 



+ 2m(0i - 02)^(a - 4M2) 
')i = tan^^(u7) 02 = ta,n^^{w/v) 



(116) 



/C > and v = w (0 < a < AM^) 



l3ia\b,c) = ^ 2Cl2(20i)- 012(201- 202) - 012(201 + 202) 



+ 202 log 



w — V 



W + V 



^1 = tan"^(l/tf;) 



2i7r0i6'(|w| — v) 
02 = tan~"^(l/iJ) 



(117) 



Note that the prefactor of i in the above equations cancels against the i from p = ip 
in the prefactor of Eq. ljllUl . so that the terms involving the Olausen function Ol2(x) 
(discussed in Appendix C) contribute to the real part of Jm- 



Case JC <0 



Here x is real. 



/C < and v real <^=^ (a < or a > 4M^) 



l3{a\b, c) = Re (^112(2:1) - Li2(zi) + 112(^2) - Li2(2;2)) 
+ 2in((pia{a)9{\x\-v) + (p29{a-4:l\/P 



^1 



^log 



X + 1 



X — 1 



V?2 



^log 



(x + f )(x — 1) 



(x — f )(a; + 1) 



(118) 



/C < and w = w (0 < a < 4M2) 



/3(a|&,c) = 2Re (112(21) - Li2(^r)) 

X + 1 X — 1 



X — IV 



X + iv 



(119) 
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Several features of these results deserve comment. 

First, in the JC > 0,v = iv case, there are anomalous thresholds which give rise 
to a nonzero imaginary part and a diverging real part. As seen in Eq. ()117|l . these 
anomalous thresholds occur in l3{a\b,c) (and similarly for J3(6|c, a) and Is{c\a,b) by 
permutation) when 

\w\ = V abc = AM^JC. (120) 

There will be a nonzero imaginary part for > U <^==^ abc > 4M^/C. Note that since 
here 4M^ > a > and /C > 0, we must have b,c> 0. Let us now look at some special 
cases: 

• a = b = c Here the condition for an anomalous threshold reduces to a > 3M^, 
which was found in 



• b = c This leads to (a/A/P) = {b/M^){4: - (b/M^)), which is possible only if 
b < 4M2. 

There are also anomalous thresholds for the case of /C < and a > 4M^. For example, 
for the mixed signature symmetric case a = b = —c > 0, there is a discontinuity in 
the real part of Jm{0'i a? and a divergence in the imaginary part at a = 5M^, as 
seen in Fig.(jn)). Anomalous thresholds were analyzed long ago |^ |49j . 

In the case /C > 0,t; real, there is an imaginary part above threshold, a > 4M^, 
which vanishes in the massless limit — oo. 

In [32], the authors find an interesting geometrical interpretation and derivation 
of the triangle integral (and higher n-point integrals). 

In the symmetric limit a = b = c, the above results reduce to those given in 
Eqs. (55-62) of [22] • 

In the massless limit, we obtain 



/C > 



J(a,6,c) = --(ci2(2(/.„) + Cl2(20^) + Cl2(20^)) 



/C < 



P 

arctanf — ), etc. (1211 
a 



Jia,b,c) = -^^Clh2(20„) + Clh2(20^)+Clh2(20g 

+ m(j)aO{a) + m(f)fj9{b) + iinp^O^c 
a + p 



- log 

2 ^ 



a — p 



etc. (122) 
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where 

Clh2(20,) for ab>0 \ 
Clh2(20„) - ^ Aah2(20„) for ab<0 J ' ^^^^^ 

and similarly for Clli2(2</)^) when {be > 0,bc < 0) and Clh2(2</)^) when (ca > 
0,ca < 0). 

The results for the massless case are well known [U] [31] [33] jSH], although the notation 
is non-standard. Here we have adopted the notation of [3H] by using the hyperbolic 
Clausen function Clh2(x), and alternating hyperbolic Clausen function AClh2(x), 
which are discussed below. 



Appendix C: Special Functions 

Here we collect some useful results, mainly taken from [SU]. The dilogarithm 
function is defined for complex z by 

U2{z) = - r dx log iln^. (124) 

In order to find the real and imaginary parts of this function, one should first ensure 
that the modulus is less than unity by judiciously using 

Lhiz) = -Li2(l/z) -^-l log' i-z). (125) 

The notation Li2(r, 9), with two arguments, is used for the real part of Li2(re''^). For 
modulus less than unity, r < 1, we have the integral representation 

riog-(l-2xcos^ + x^) 
2 Jo ^ 

The imaginary part for r < 1 is 

Im (Li2(re'')) = T logr + ^ (Cl2(2^) + Cl2(2T) - Cl2(2^ + 2T)) 

^ T / r sin 6 \ , ^ 

T = tan^M . (127 

\1 — r cos 9 J 

In particular, 

Im (Li2(e'')) = Cl2(^) and Cl2(^) = ^ (Li2(e'^) - Li2(e-'^)) . (128) 
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The Clausen function frequently appears in the triangle integral and has the following 
representations : 



^, / N , 1^ . Vi v~^sinnx 

Cl2(x) = - / dylog|2sin|| 

-^0 n=l 



Furthermore, Cl2(x) is odd, Cl2(— x) = — Cl2(x), satisfies periodicity, Cl2(a; + 2n7r) = 
Cl2(a;), and a duplication formula Cl2(2a;) = 2Cl2(x) + 2Cl2(x — vr). Many other 
properties can be found in [201 the some are conveniently summarized in the 
appendix of 

We have used the notation of Luj2Hl, used the hyperbolic Clausen function, 
Clh2(x), and alternating hyperbolic Clausen function, AClh2(a), defined by the inte- 
gral representations 

r y 

Clh2(a;) = — (i|/log |2sinh-| 
Jo 2 

AClh2(x) = - / rfylog|2cosh-|. (129) 
Jo 2 

These can also be written as 

Clh2(a:) = ^Re(Li2(e^)-Li2(e-^)) 

AClh2(a;) = ^Re(Li2(-e") -Li2(-e-")) (130) 

in analogy with Eq. p28|l . 

Finally, some elementary relations which are used often include (for x, y real) 

arg (x + iy) = ta,n~^ ^— j + iT6{—x)a{y) 

(X ~\~ if \ 
I + 7ra(x)9(xy — 1) 
1 - xyj 

tan~^(x) +tan~^ (1/x) = (t(x)|. (131) 

where cr{x) = x/\x\ is the sign function and the step function 6{x) = (cr{x) + l)/2 
should not be confused with the angle 6. 



Appendix D: Form Factors with Dimensional Re- 
duction (DRED) Regular izat ion 

Here we discuss the form factors regularized using dimensional reduction (DRED) 
in integer number of dimensions d^, defined analogously to the usual = 4 DRED 
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scheme. This could be used for cIr = Q or dji = 10 theories, but of course we mainly 
have in mind the four-dimensional case. 

It is easy to see that the quark and scalar contributions are unchanged from 
DREG, and only the gluon contribution is different. This is most easily expressed in 
terms of the modified sum rule 

FGiDRED)+AFQ + (10 -dR)Fs = 0, (132) 

which implies Fg{DRED) = Fg{DREG) + (dn - d)Fs. Expanding d = dn - 2e 
around the real number of dimensions dji leads to Fg{DRED) = Fg{DREG) + 
2eFsi which makes manifest the role of the 2e adjoint DRED "ghosts" which preserve 
supersymmetry. 

In four dimensions we have 

Fg{DRED) + AFq + 6F5 = 0. (133) 

Since only the A form factors have a UV divergence in four dimensions, only these 
form factors will be changed when using DRED: 

Wd(^i2(G)) = 5dred(Ai2{G)) = (134) 

In the symmetrized physical ± basis we have 5dred{A.q{G)) = —^j^, and all other 
form factors are unchanged. 



Appendix E: Quark and Squark Mass Corrections 

Here the corrections to the form factors due to fermion and scalar masses will 
be given. The massive quark (MQ) contributions were first obtained in (HHl, and we 
obtain exactly the same results. To our knowledge, the squark contributions, either 
massless or massive (MS), have not yet appeared in the literature. 

First, the well known formulas for the scalar and fermion self-energies are repro- 
duced in our notation : 

2JiM-{d-2)TM \ 
a{l -d) J 

with the integrals Jim, Tm given in Appendix A. These yield the massive fermion and 
scalar contributions to the longitudinal form factors through Eq. ()38|) . 



ni(MQ) 



d-2 
1 - d' 



Jim + 2M' 
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Notice that the terms proportional to in the above two equations are equal 
up to a factor of —2. After explicit calculation, it was discovered that the scalar mass 
correction terms {6ms) are just minus one-half of the quark mass correction terms 
(Smq) for all form factors. Thus, for generic form factor F 



F{MS) 
F{MQ) 



FiS) +6MsiF) 

M 

F{Q) +SMQiF) 

M 



(136) 



The notation F{S) simply means to take the appropriate massless result for the 

M 

form factor F, as given in section 3, and replace J — > Jm, Ji Jim, J2 -hM, 
J3 JsM everywhere. 

Because of the relation 6mq = —'^^ms, we need only write either the fermion 
or scalar mass correction terms explicitly. Here we choose the scalar contributions, 
which for the transverse form factors are 



Sms{F 



12 



2M2 



/C-3/?7 
d-2 



V - 2a/C - 72(3a + 2-f - /3) 



M 



a(/3-7)(c/-l) 



Jim (137) 



V - 2a/C - /52(3a + 2/5-7) 
h{(3-^){d-l) 



Jim — 



2c 



■J 



3M 



(2-rf)/Ca, 
2{d-l)ab' 



■M 



and 



Sms{H) 



3V 
d-2 



J, 



M 



/C-2a7 
d- 1 



J 



IM 



K-2aB K. 
—J2M - 



d 



d 



2/37 . ^ lC{d 



2{d 



-T, 



M 



(.138) 



The results in the physical ± basis can be obtained from those of the LT basis 
through the use of Eq. lfTBj) . and are included here for completeness. 



M2 

Sms^Au) = — 



ca 7 Jim + /3J2m + cJsm 
:Jm H 



d-2 



d 



(139) 



hisiB 



12; 



M2 

Ic 



[a-h)a {2a + ^)JiM - {2a + P)J2M + {P - l)J3M 
-•Jm H ; ; 



d-2 



(140) 
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his{C_ 



12 



W + a^c + 3/3^(a - 7) + 37^(a - (3) 
d-2 



J. 



M 



+ — — 7^ Jim 



a{d- 1) 

+ a\-f + p\5a + 77 + 4/j) 
bid-1) 



J2M 



4c2 + 2/?7 - /C ^ (rf - 2)/C(a2 + 2/^7 + 3/C) ^ 
+ ; : J'iM H Ti 7T7r~i -'A/ 



d-l 



{d-l)2ab 



(141) 



^A/5(^12) 



M2 
4/C2 



(/C + 3a2)(Q _ 5) /c + 4^2 ^ 2a7 
— Jm + - 



d-2 



d 



IM 



(142) 



]C + Aa'^ + 2aP (2/57 - 3/C) (a - 6) (d-2)/C(a-6) 
-<>'2M H — — 7\ JsM H :rr, — ~ 



d 



c{d 



2{d-l)c 



5ms{H) 



M2 
4)C 



2V + ahc 



2-d 



J 



M 



(rf-2)/C 
2(d- 1) 



T. 



M 



/C-2a7 
d- 1 



lA/ 



/C - 2a/5 /C - 2/57 

--'2M H ; ; — 



rf- 1 



d-l 



(143) 



(^Afs(>S') 



M2 



4/C2 



3(a — h)(h — c)(c — a) 
2 - 



J, 



A/ 



(4a2 + 2a7-3/C)(fe-c) 
a{d-l) 



IM 



(452 + 2a/3 - 3/C) (c - a) (4c2 + 2/37 - 3/C) (a - 6) , 

+ 7T-, 7^ J2M H Tl 7\ "^SAf 



6(c/- 1) 

{d - 2){a - h){h - c){c - a)lC ^ 
2{d- l)ahc " 



c{d-l] 



■M 



(144) 



It is straightforward to see that all of the correction terms are ultraviolet finite. 
The relation —25ms = ^mq is necessary for the preservation of the form of a 
quark/scalar sum S^q = 2Fs + -Fq, which is equal to -j^T^qg using the results of 



2Fs\,^ + Fc 



lAr- 



so that this quantity has 



section 3. In other words 2Fms + Fmq - 
no correction terms proportional to M^. 

However, the relations between massive gauge bosons and massive fermions and/or 
scalars will be different, since the gauge bosons eat a degree of freedom to acquire 
mass. Consider the contribution of a massive gauge boson to the gauge-invariant 
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gluon self-energy 



Ui{MG) 



J 



IM 



]d — dR — 7 
2{d-l) 



2M2 {dn 



a {d — 1) 



JiM + ^{2-d)TM 



(145) 



where as before da = d in dimensional regularization (DREG) and dn = 4 (or the 
real integer number of dimensions) in dimensional reduction (DRED). From this and 
Eq. p35|) we deduce 



Ui(MG) + 4ni (MQ) + (9 - dR)Ui{MS) = 
and thus the massive = 4 sum rule becomes 



(146) 



Fmg + ^Fmq + (9 - dR)FMs = (147) 

for the longitudinal form factors. It can also be shown that this holds for the trans- 
verse form factors and so the results of this paper also give the contributions of massive 
internal gauge bosons. Proving this involves detailed analysis of the vertices and dia- 
grams that contribute to the triple gluon vertex when the PT/BFMFG is applied to 
a spontaneously broken gauge theory that leaves a non-abelian subgroup (of gluons) 
intact. This can be done following a pinch-technique route similar to ^SICH]- Due to 
the equivalence of the PT and BFMFG, it is more convenient to follow the BFMFG 
route similar to [5T]. For example, in SU{5) GUTs the colored superheavy X and 
Y gauge bosons give a contribution which satisfies the above massive sum rule. We 
should emphasize that these sum rules are simply a convenient way of relating the 
contributions of various spin particles, all with hypothetical mass M, but no assump- 
tion is made about the actual masses for a given theory under consideration; the sum 
rules are entirely stripped of color factors. 



^This is also the contribution of to the photon self-energy. 
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